THE CONTRACTION RATE IN THOMPSON METRIC OF 
ORDER-PRESERVING FLOWS ON A CONE - APPLICATION TO 
GENERALIZED RICCATI EQUATIONS 



STEPHANE GAUBERT AND ZHENG QU 

Abstract. We give a formula for the Lipschitz constant in Thompson's part metric 
of any order-preserving flow on the interior of a (possibly infinite dimensional) closed 
convex pointed cone. This provides an explicit form of a characterization of Nussbaum 
concerning non order-preserving flows. As an application of this formula, we show 
that the flow of the generalized Riccati equation arising in stochastic linear quadratic 
control is a local contraction on the cone of positive definite matrices and characterize its 
Lipschitz constant by a matrix inequality. We also show that the same flow is no longer a 
contraction in other natural Finsler metrics on this cone, including the standard invariant 
Riemannian metric. This is motivated by a series of contraction properties concerning 
the standard Riccati equation, established by Bougerol, Liverani, Wojtowski, Lawson, 
Lee and Lim: we show that some of these properties do, and that some other do not, 
carry over to the generalized Riccati equation. 



I. Introduction 

The standard discrete or differential Riccati equation arising in linear-quadratic control 
or optimal filtering problems has remarkable properties. In particular, Bougerol |Bou93j 
proved that the standard discrete Riccati operator is non-expansive in the invariant Rie- 
mannian metric on the set of positive definite matrices, and that it is a strict contraction 
under controllability /observability conditions. Liverani and Wojtowski |LW94] proved 
that analogous contraction properties hold with respect to Thompson's part metric. 
These results, which were obtained from algebraic properties of the linear symplectic 
semigroup associated to a Riccati equation, are reminiscent of Birkhoff's theorem in 
Perron-Frobenius theory (on the contraction of positive linear operators sending a cone 
to its interior |Bir57j ). Lawson and Lim |LL07j generalized these results to the infinite 
dimensional setting, and derived analogous contraction properties for the flow of the 
differential Riccati equation 

(1) P = A'P + PA- PEP + Q, P(0) = G , 

where A is a square matrix, S, Q are positive semi-definite matrices, and G is a positive 
definite matrix. Moreover, Lee and Lim |LL08j showed that the same contraction prop- 
erties hold more generally for a family of Finsler metrics invariant under the action of 
the linear group (the latter metrics arise from symmetric gauge functions). 
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It is natural to ask whether contraction properties remain valid for more general equa- 
tions, like the following constrained differential Riccati equation, 

P = A'P + PA + C'PC + Q 
, s -{PB + C'PD + L')(R + D'PD)-\B'P + D' PC + L), 

[ ] P(0) = G 

R + D'PD positive definite, 

which has received a considerable attention in stochastic linear quadratic optimal control. 
The equation ([2]) is known as the generalized Riccati differential equation (GRDE) or as 
the stochastic Riccati differential equation. Up to a reversal of time, it is a special case of 
the Backward stochastic Riccati differential equation, which have extensively studied, see 
in particular |YZ99l CLZ98, RCMZOlaJ. The reader is referred specially to the monograph 
by Yong and Zhou [YZ99] for an introduction. Even for the simpler Riccati equation ([1]), 
contractions properties have not been established when the matrices Q, E are not positive 
semi-definite, whereas this situation does occur in applications [McE07j. 

In this paper, motivated by the analysis of the generalized Riccati equation, we study 
the general question of computing the contraction rate in Thompson's metric of an ar- 
bitrary order-preserving (time-dependent) flow defined on a subset of the interior of a 
closed convex and pointed cone in a possibly infinite dimensional Banach space. Recall 
that the order associated with such a cone C is defined by x^y-^y — x&C, and that 
the Thompson metric can be defined on the interior of C by the formula 

d T (x,y) := log(max{M(x/y), M(y/x)}) 

where 

Mix/y) := inf{t G R : ty ^ x} = sup ^4 , 

^ec* ip(y) 

and C* denotes the dual cone of C. More background can be found in §2.11 
Our first main result can be stated as follows. 

Theorem 1.1. Assume that the flow of the differential equation xit) = <f)(t,x(t)) is order 
preserving with respect to the cone C, and let U denote an open domain included in the 
interior of this cone such that XU C U holds for all A G (0, 1]. Then, the contraction rate 
of the flow over a time interval J, on the domain U, with respect to Thompson metric, 
is given by the formula 

(3) a := — sup M((_D0 s (x)x — <p(s, x))/x) . 

se J, x&a 

Here, D(f) s (x) denotes the derivative of the map (s,x) i— > <f)(s,x) with respect to the 
variable x. We make some basic technical assumptions (continuity, Lipschitz character 
on the function with respect to the second variable) to make sure that the flow is 
well defined. We refer the reader to Section [3] for more information, and in particular to 
Theorem 13.51 below, where the definition of the contraction rate can be found. 

The idea of the proof is to construct a special flow- invariant set, appealing to a gener- 
alization due to Martin |Mar73] of theorems of Bony |Bon69] and Brezis |Bre70] on the 
geometric characterization of flow invariance. Formula 03) should be compared with re- 
sults of Nussbaum, who studied the more general question of computing the contraction 
rate of a non- necessarily order-preserving flow in Thompson metric |Nus94j , and obtained 
an explicit formula in the specical case of the standard positive cone. However, this for- 
mula, valid for non order-preserving flows, appears to have no natural generalization to 
abstract cones (although some reasonably explicit conditions can be given in the special 
case of symmetric cones, we leave this for a further work, see also Section [6] below for 
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a special case). In addition, Nussbaum's approach, which relies on the Finsler structure 
of the Thompson metric, is widely applicable in its spirit but leads to different technical 
assumptions, including geodesic convexity. See §4.11 for a detailed comparison. 

Then, we show, in Section I4.2[ that the contraction results of Liverani and Woj- 
towski |LW94] and of Lawson and Lim |LL07] concerning the standard Riccati equa- 
tion (P) with positive semi-definite matrices £, Q, as well as new contraction results in 
the case when E is not positive semi-definite, can be recovered, or obtained, by an appli- 
cation of Formula fl3]). This provides an alternative to the earlier approaches, which relied 
on the theory of symplectic semigroups. This will allow us to handle as well situations 
in which the symplectic structure is missing, as it is the case of the generalized Riccati 
differential equation. 

Our second main result shows that the flow of the generalized Riccati differential equa- 
tion is a local contraction in Thompson metric. 

Theorem 1.2. Assume that the coefficients of the generalized Riccati differential equa- 
tion 02]) are constant, and that the matrix ( ^ L R ) is positive definite. Then, the flow of 
this equation is a strict contraction on the interior of the cone of positive definite ma- 
trices, and this contraction is uniform on any subset that is bounded from above in the 
Loewner order. 

This theorem follows from Theorem 15.31 in Section where an explicit bound for the 
contraction rate on an interval in the Loewner order is given. We shall also see in Sec- 
tion |5]that the flow of the generalized Riccati equation is no longer a uniform contraction 
on interior of the cone, which reveals a fundamental discrepancy with the case of the 
standard Riccati equation. Then, motivated by earlier results of Chen, Moore, Rami, 
and Zhou (see |RC MZ01b| and |RZOO] ) on the asymptotic behavior of the GRDE, we 
identify (Theorem 15 .6 j) different assumptions under which a trajectory of the GRDE con- 
verges exponentially to a stable solution of the associated Generalized Algebraic Riccati 
Equation (GARE). We also establish (Section l5.4p analogous results concerning the dis- 
crete time case. Then, we give a necessary and sufficient condition (Proposition 15. lOj) for 
the generalized discrete Riccati operator to be a strict global contraction. 

Finally, in Section [6j we establish the following negative result, which shows that the 
Thompson metric is essentially the only invariant Finsler metric in which the flow of the 
GRDE is non-expansive for all admissible values of the matrix data. 

Theorem 1.3. The flow of the generalized Riccati differential equation is non-expansive 
in the invariant Finsler metric arising from a symmetric gauge function, regardless of the 
parameters (A, B, C, D, L, Q, R), if and only if this symmetric gauge function is a scalar 
multiple of the sup-norm. 

In particular, the flow of the GRDE is not non-expansive in the invariant Riemannian 
metric, showing that Bougerol's theorem on the contraction of the standard discrete 
Riccati equation does not carry over to the GRDE. 

2. Preliminaries 

2.1. Thompson's part metric. We first recall the definition and basic properties of 
Thompson's part metric. 

Throughout the paper, X is a real or complex Banach space with norm | • | . Let X* be 
the dual space of X . For any x G X and q G X*, denote by {q, x) the real part of q(x): 

(q, x) = Re q(x) . 
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Let C C X be a closed pointed convex cone, i.e., aC C C for a G M + , C + C C C and 
C R (— C) = 0. The dual cone of C is defined by 

C* = { z G : (z,x) ^ Vx G C} . 

We denote by Co the interior of C. We define the partial order ^ induced by C on X by 

x^y^y — xEC 

so that 

x ^ y ^ (z,x) ^ (z, y) , V2 G C*. 
We also define the relation <C by 

x <^ y y — x E Co . 

For x ^ y we define the order intervals: 

[x,y] := {z G A'jx ^ z < y}, (x,y) := {2 G ATjrc < z <?/}. 

For a; G <Y and y £ Co, following [Nus88] . we define 

, M(x/y) := mf{t eR:x^ty} 
^ ' m(x/y) := sup{t Gl:i) ty} 

Observe that since y G C , and since C is closed and pointed, the two sets in (TJ]) are 
non-empty, closed, and bounded from below and from above, respectively. In particular, 
m and M take finite values. 

Definition 2.1. The Thompson part metric between two elements x and y of Co is 

(5) d T (x,y) :=\og(max{M(x/y),M(y/x)}) . 

It can be verified that dx(-, ■) defines a metric on Co, namely for any x,y,z G C we 
have 

d T (x,y)^0, d T (x,y) = d T (y,x), d T {x, z) ^ d T {x,y) + d T {y, z), d T {x, y) = <^> x = y. 

A sufficient condition for Cq to be complete with respect to dr{-,-) is that C is a normal 
cone, see |Tho63j . We shall consider specially the case in which C is the cone of n x n 
positive semi-definite matrices. Then, it can be checked that, for all A,B G Co, 

M(A/B) = max logAj, , 

where Ai, . . . , A„ are the eigenvalues of the matrix B~ X A (the latter eigenvalues are real 
and positive) so that the Thompson metric d? can be explicitly computed from (JSJ). 

2.2. Characterization of flow invariant sets. We next recall some known results on 
the characterization of flow-invariant sets in terms of tangent cones, which will be used 
to characterize order-preserving non-expansive flows in Thompson's metric. 

In the sequel, J = [0, T) C K is a possibly unbounded interval, V C X is an open set, 
and 4>(t, x) is a function from J x T> to X . For x G X and S C X we define the distance 
function: 

d(x, S) = inf{|x — y\ : y G S}. 
We study the following Cauchy problem: 

x(t) = <p(t,x(t)), 

x(s) = Xq. 



(6) 



By a solution of ([6]) on [s, a) C J we mean a continuously differentiable function t 1— > 
x(t) : [s, a) — > V such that x(s) = Xo and x(t) = <fi(t, x(t)) for all t G [s, a). 
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Let S be a closed subset of X. We say that the system [S fl D, 0) is flow-invariant if 
every solution of (E]) leaves S invariant, in the sense that for any s G J and Xq G S fl T>, 
the solution x(t) must be in G 5 fl V, for all t G [s, a). 

Characterizations of flow invariant sets go back to the works of Bony |Bon69] and 
Brezis |Bre70] . Several improvements, together with extensions to the infinite dimensional 
case can be found in |Red72] . |Mar73] . |Cla75] and |RW75] . We shall actually need here 
an immediate consequence of a theorem of Martin |Mar73] . 

Theorem 2.2 (Theorem 1 of [Mar73j). Suppose that the following conditions hold: 
(CI) is a continuous function on J x T>; 

(C2) For every closed bounded set K dV, there is a constant L > such that 

\<f>{t,x) -<f>{t,y)\ ^ L\x-y\, Vt G J,x,y G K ; 

(C3) For allt G J and x G S DV, 

, d(x + h6(t,x),S r\V) 

hm = . 

hiO h 

(C4) S is convex. 
Then the system (S (1 V, <ft) is flow-invariant. 

It is not difficult to prove that for i 6 5nP, v E X and sufficiently small h > 0, 

d(x + hv, S) = d(x + hv,SnV). 

Thus, Condition (C3) is equivalent to: 
(C5) For allt G J and x G <S n V, 

,. d(x + h(j>(t,x),S) 

hm = 0. 

hiO h 

Condition (C2) is a local Lipschitz condition for the function 0, with respect to the 
second variable. Condition (C3) is a tangency condition (the vector field should not 
point outward the set S fl V). 

Definition 2.3 (Tangent cone |Cla75j ). The tangent cone to a closed set S C X at a 
point x G S, written Ts(x), is the set of vectors v such that: 

(7) hminf^±^ = 0. 

Ho h 

Remark 2.4. When S is a closed convex, we know that the limit in (j7j) exists. Thus, 
Condition (C5), equivalent to (C3) in Theorem 12.21 can be replaced by: 

(C6) For allt G J and x G S n V, 

<f>(t,x) G T s {x). 

Besides, this definition coincides with the one in convex analysis, i.e., 

Proposition 2.5 (Proposition 5.5, Exercice 7.2 |CLSW98j). Let S be a closed convex set 
of X , then, 

T s (x) = c\{v : 3A > with x + Xv G S}, Vx G S. 

Now the flow-invariance can be checked by verifying if 0(t, x) lies in the tangent cone 
of S. To this end, we need to compute the tangent cone at each point of S. In some 
cases the tangent cone can be expressed in a simple way: 
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Lemma 2.6 (Exercise 2.5.3 [CLSW98J). LetSi,S2 C X be closed subsets, X — (^Xi^X 2 ^j G 
S\ x S 2 . Then 

T Sl xs 2 (x) = T Sl (xi) x T S2 (x 2 ). 

We shall consider specially S = C. Then, using Proposition 12.51 and the Hahn-Banach 
separation theorem, one can show that 

(8) T c {x) = {v\(q,v) ^ if q G C* and (q,x) =0}, x G C. 

3. The contraction rate in Thompson metric of order-preserving flows 

3.1. Preliminary results. From now until the end, the function <f)(t,x) is assumed to 
be continuous on J x D and Frechet different iable to x. The derivative of <ft with respect 
to the second variable at point (t,x) is denoted by D<f) t (x). We also assume that the 
derivative is bounded on any closed bounded set, i.e., for any bounded set K C V, there 
is a constant L such that: 

\D(j) t (x)\ ^ L, MteJ.xeK. 

Therefore Condition (CI) and Condition (C2) are both satisfied. The existence and 
uniqueness of the solution of ([6]) follow from the Cauchy-Lipschitz Theorem. We then 
define the flow M.'(-) associated to the system by: 

Ml(x Q ) = x(t), s G J, t G [s, a) 

where x(t) : t G [s, a) is the maximal solution of ([6]). (Note that in general the flow is 
defined only on a subset of J x J x V.) For each open subset U C V and initial value 
x G U we define tu(s, x ) as the first time when the trajectory leaves U, i.e., 

tu(s,x ) = sup{6 G (s,a)\Ml(x ) G U, Wt G 

When is independent of time t, we denote simply 

M t (x ) := M*(x ), tu(x ) := t w (0,x ), Va;o EV,UcV. 
By uniqueness of the solution, the flow has the group property: 

Ml(x ) = M^HM^xq)), V ^ s ^ h <: t < t v {s,x ) . 

Definition 3.1 (Order-preserving flow). Let U be an open subset of T>. The flow M'(-) 
is said to be order-preserving on U if for all Xi,x 2 6W such that Xi ^ x 2 , 

M l s (xi) ^ M*(x 2 ), V0^s^(< t u (s, xi) A t w (s, x 2 ). 

Definition 3.2 (Non-expansiveness and contraction). Suppose that C C T>. The flow 
M.'(-) is said to be contractive on C with rate a > in Thompson metric if for all 
Xi, x 2 G Co, 

^(M^xxXM^)) ^ e~ a{t - s) d T (xi,x 2 ), V0O^< *cb(s,a;i) A t Co (s, x 2 ) 
If the latter inequality holds with a = 0, the flow is said to be n on- expansive. 

In the following, our primary goal is to characterize non-expansive order-preserving 
flows in Thompson part metric. We shall need the following proposition, which provides 
a characterization of monotonicity in terms of the function <fi. The equivalence of the first 
two assertions was proved in |RW75j . 

Proposition 3.3 (Compare with Theorem 3 in |RW75j ) . LetU be an open subset ofD. 
The following conditions are equivalent: 

(a) The flow M.'(-) is order-preserving on U. 

(b) For all s G J and Xi, x 2 G U such that x\ ^ x 2 , <f>(s, Xi) — (f)(s, x 2 ) G Tc{x% — x 2 ). 
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If 14 is convex, then the above conditions are equivalent to: 
(c) For all s E J x El4 and v G C, 

(9) (g, D4> s {x)v) > 0, Vg G {g G C* : (g, u) = 0}. 

Proof. We only need to prove the equivalence between (b) and (c), since the equivalence 
between (a) and (b) follows from [RW75] . In view of flSJ), Condition (b) is equivalent to 
the following: 

for all s G J and Xi,x 2 ElA such that X\ ^ x 2 , 

(q, (f)(s, xt) - (f)(s, x 2 )) ^ 0, Vg G {q G C* : (g, Xi - x 2 ) = 0}. 

Now suppose that (b) is true. Then for any s G J, x G 14 and any v G C, there is 5 > 
such that for any ^ e ^ 5 

(g, 0(s, x + eu) - (f){s, x)) ^ 0, Vg G {g G C* : (g, v) = 0}. 

Since <p is differentiable at point x, dividing by e the latter inequality, and letting e tend 
to 0, we get 

(g, D<j) s {x)v) > 0, Vg G {g G C* : (g, u> = 0}. 

Next suppose that Condition (c) holds. Fix any s G J and Xi, x 2 ElA such that Xi ^ x 2 . 
Fix any g G C* such that (g, Xi — x 2 ) = 0. Define the function g : [0, 1] — >• M by: 

g(A) = (g, 0(s, Axi + (1 - A)x 2 ) - 0(s, x 2 )). 

Then we have g(0) =0 and in view of convexity of 14 and ([91 , 

g'(A) = (g, D0 s (Ax! + (1 - A)x 2 )(x! - x 2 ))) ^ 0, V0 < A < 1. 

A standard argument establishes that: 

(10) g{l) = {q,(j>{s,x 1 )-<f>{s,x 2 ))>0. 

Since s,xi,x 2 and g are arbitrary, we deduce Condition (b). □ 

3.2. Characterization of the contraction rate in terms of flow invariant sets. 

The following is a key technical result in the characterization of the contraction rate of 
the flow. 

Proposition 3.4. Let 14 CD be an open set such that X14 C U for all A G (0, 1]. If the 
flow M.'(-) is order-preserving on 14, then the following conditions are equivalent: 

(a) For all x EU and A ^ 1 such that Xx G 14, 

Ml(Xx) < X e ~ a{t ~ a) M l s (x), < s < t <t u (s,x) A t M (s, Ax). 

(b) For all s G J and x Ell, 

D(p s (x)x — 4>(s,x) ^ —ax. 

(c) For a// x,y El4 and X ^ 1 sucA t/iat y ^ Ax ; 

Mliy) ^ A e " a(t " s) M s *(x),0 ^ s < t < t u (s,x)At u (s,y). 

Proof. Suppose Condition (a) holds. Let x be any point in U. Fix any A > 1 such that 
Ax EU, we must have: 

(11) Ml(Xx) < A e " a(t_s) M*(x), < t < t u (s, x) A tu(s, Xx). 
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where it must be the case that tu(s,x) A t u (s,Xx) > 0. Since the terms on both sides 
of (II ip coincide when t = s, taking the derivative of each of these terms at t = s, we 
obtain 

(12) <f)(s,Xx) ^ X<p{s,x) - a{XlnX)x 

Since this inequality holds for all A ^ 1 such that Xx G U, with equality for A = 1, the 
derivation of the two sides of the above inequality at A = 1 leads to: 

Dcf) s (x)x — (f>(s, x) ^ —ax 

for all x G U. Condition (b) is deduced. 

Now suppose that Condition (b) is true. We shall derive Condition (c) by constructing 
an invariant set. Denote: 

X := X x X x R, 
V:=UxUx M+\{0}, 
S := {(xi, X2, A) G X : x 2 ^ Xx\, A ^ 1}. 
Define the differential equation on T>: 

( x x \ / <f>(t, Xl ) \ 

(13) x 2 = ${t,x 1 ,x 2 ,X) := 4>(t,x 2 ) 

\ A / V -aXlnX J 

It is not difficult to see that Condition (c) is equivalent to the flow-invariance of the 
system (S flD,!). It would be natural to show directly the latter flow-invariance by 
appealing to Theorem \2.2\ but the set S is not convex, making it harder to check the 
assumptions of this theorem. Therefore, we make a change of variable to replace S by a 
convex set. 

Define the smooth function F : X — > X by: 

F(x\, x 2 , A) = (xi, Axi — x 2 , X — 1), V(xi, x 2 , A) G X. 

Denote 

S' = X x C x R+. 

By Lemma [276]. for (yi,y 2 , k) G S', 

T S '(yi,y2,K) = X x T c (y 2 ) x T r +(k). 

Observe that S = {x G X\F(x) G S'} and that F has a smooth inverse G : X — > X given 
by: 

G(y h y 2 , k) = (y h {k + l)yi -y 2 ,K + l). 

Therefore F(V) = G^iV) is an open set. Let (yi,y 2 , k) = F(xi,x 2 , A) and consider the 
system: 

(yi,y2,K)' = ^(t,y 1} y 2} K) 

where 

[ <t>(t, yi) \ 
${t, y u y 2 , k) := -a(K + 1) ln(/t + l)y x + (k + l)0(t, yx) - 0(t, (k + - y 2 ) 

\ -a(K + 1) ln(K + 1) / 

One can verify that the invariance of the system (iSnP, $) is equivalent to the invariance 
of the system (S' n ^). 

Now the function $ : /x -F('E') — >■ ^ defined as above is continuous and differentiable 
to the second variable with bounded derivative on bounded set. Besides S' is convex. 
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By applying Theorem \2.2\ the system (<S' fl F(t>),^) is flow-invariant if the following 
condition is satisfied: 

(14) y(s,y u y 2 ,K) eT s ,( yi ,y 2 ,K), Vs G J, ( yi ,y 2 , k) G <S' D F{V). 
That is, for any y 2 , K ) G 5' fl .F(X>) and s E J, 

0(s,yi) G # 

-ar(« + 1) ln(« + 1)^ + (k + l)0(s, j/i) - (k + l)yi - y 2 ) G T c (y 2 ) 
-a(« + 1) ln(re + 1) G T r +(k) 

It suffices to check the second condition because the others hold trivially. By applying 
the bijection F, this condition becomes: for any s G J and (xi,x 2 , A) G S fl P, 

— aA In Axi + A0(s, xi) — 0(s, x 2 ) G Tc(Xxi — x 2 ). 

Let any s G J, xi,x 2 G U and A ^ 1 such that x 2 ^ Axi. Let any q G C* such that 
(g, Axi — x 2 ) = 0. By dHJ) we only need to prove: 

(15) (g, — aA In Axi + X(f>(s, Xi) — (f>(s, x 2 )) ^ 0. 

By the assumptions, we know that A _1 x 2 G U. Then, it suffices to prove: for any 
Xi,x 2 G U such that X\ ^ x 2 , let q G C* such that (g,Xi — x 2 ) = 0, then for any A ^ 1 
such that Ax 2 G W we have: 

(q, — aA In Axi + A0(s, Xi) — 0(s, Ax 2 )) ^ 0. 

Define the function / : [1, A] — >■ R by: 

/( r ) = (<?> -alnrxi + 0(s,Xi) - r~ l (j)(s,TX 2 )) 

Notice that the function / is well defined on [1,A]. By hypothesis of monotonicity and 
Proposition 13.31 

/(l) = (g,0( S ,x 1 )-0( S ,x 2 ))^O. 
Differentiating / gives, for all r G [1, A], 

f(r) = (q, -r~ 1 ax 1 + t~ 2 0(s,tx 2 ) - t' 1 D<P s (tx 2 )x 2 ) 
^ (q, —t olx\ + t~ x olx 2 ) (by Condition (b)) 
= . 

A standard argument establishes that /(A) ^ 0, and so f[T5"j) is proved, whence the flow- 
invariance of [S flP,$), which is exactly Condition (c). Finally, Condition (a) follows 
from Condition (c) by considering y = Ax. □ 

We next state the main results. Recall that J = [0,T) C R. 

Theorem 3.5 (Contraction rate). Assume that <fi is defined on J x U where U <Z Cq is 
an open set in the interior of the cone such that \U C U for all A G (0, 1]. If the flow 
M'(-) is order-preserving onU, then the best constant a such that 

(16) d T (Ml(xi), M*(x 2 )) < e~ Q( *- s W(xi, x 2 ), «C s < t < t u (s, xi) A t u (s, x 2 ) 
holds for all x\, x 2 G U is given by 

(17) a := — sup M((D(f) s (x)x — <f>(s, x))/x) . 

sG J, x£U 
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Proof. If (TT6"P holds for all X\,x 2 E U, then Condition (a) in Proposition ( 13. 4 p holds. It 
follows that the constant a must satisfy 

(18) D(f> s (x)x — <fi(s,x) ^ —ax, Ws E J, x E U. 

Now conversely if (fl8l) holds. Then Condition (c) in Proposition 13.41 holds. For any 

xi,x 2 E U, let A = e dT(xi ' X2 \ then 

M*(xi) =sC A e " Q(t " s) M*(x 2 ),0 < s<:t< t u (s,x 2 ) At u (s,xi) . 

The same is true if we exchange the roles of X\ and x 2 , and so, ffl~6]) holds for all x±, x 2 E U. 
Consequently the best constant a such that (TIB]) holds for all Xi,x 2 E U must be the 
greatest constant a such that ffTgj) holds, which is precisely ffTTp . □ 

Now we get a direct corollary. 

Theorem 3.6. Suppose that <fi is defined on J x Cq. Let a G 1. If the flow is order- 
preserving on Co, then the following are equivalent: 

(a) For all x\,x 2 E Co'. 

d T (M s '(:ci),M*(x 2 )) < e- a ^d T (xi,x 2 ), ^ s ^ t < t Co (s, x x ) A t Co {s, x 2 ). 

(b) For all s E J and x E Co, 

D(f) s (x)x — 4>(s, x) ^ — ax. 

If any of these conditions holds, then the flow leaves Co invariant, i.e., for any s E J and 
x E Co, t Co (s,x) = T. 

Proof. The equivalence between (a) and (b) follows from Theorem 13.51 Now suppose 
that Condition (b) holds. Let any s E J and x\,x 2 E Co- Let t± = tc (s,x\) and 
^2 = tc {s,x 2 ). Suppose that t\ < t 2 . Then it must be the case that t\ < +oo. Thus the 
set {MJ(x 2 ) : r E [s,ti]} is compact and included in Co- Denote 

K = max{rf T (M;(x 2 ),M* 1 (x 2 ))| r E [s,ti]} < +°o 

and K = K + max{e _a ^ 1_s \ 1)^(^1, x 2 ). Note that there exists s <f <t\ such that 

d T (M:( Xl ),M*(x 2 ))>Ko, 

otherwise tc (s,xi) > t\. But for any s < r < t±, 

d T (M r s {x 1 ),M^(x 2 )) ^ d T (M:( Xl ),M r s {x 2 )) +d T (Ml(x 2 ),M*(x 2 )) 

<: e- a{ - r -^d T {x u x 2 ) + d T {M r s {x 2 ),M t ! }{x 2 )) 

The contradiction implies that t\ < t 2 is impossible. We then showed that there exists 
T E (0, +oo] such that for any s E J and x E C , tc (s, x) = T. From the group property 
of the flow action, we deduce that T — T. □ 

In the sequel we suppose that the dynamics is independent of time and study the 
convergence of an orbit of the flow to a fixed point in the interior of the cone. Let x E Co 
be such that (p{x) = 0. Let fi > 1. Denote by U the open interval (/i _1 x, fix). We look 
for the best constant a6K such that: 

(19) d T (M t (x),x) ^ e~ at d T (x, x), Va; E U, ^ t < t u (x). 
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Theorem 3.7 (Convergence rate). We assume that (ft is independent of time, defined on 
Cq and such that the flow is order-preserving on C - Let x G Cq be a zero point of (p. Then 
the best constant a such that / fT^j) holds is given by 

(20) a= inf m((-(A In X)~ l (f)(Xx)) / x) . 

Moreover, if the latter a is non-negative, then for all x G [fi^x, fix], 

(21) dr(M t (x),x) < e~ at d T (x,x), Vt ^ 0. 

Proof. Suppose that a satisfies (fT9l) . Let any A G Then Ax G W and 

M t (Ax) < A e " Qt x, < t < t w (Ax). 

Since tu(Xx) > and both sides of the former inequality coincide when t = 0, we get the 
inequality for the derivative at t = 0: 

(22) 0(Ax) ^ -aX(lnX)x , 
and so 

(23) ax ^ -(AlnA)~V(As), VI < A < p. 
Similarly, for A G (/J -1 , 1), 

A e " a 'x ^ M t (Ax), «C t < t u (Xx), 

thus 

(24) -ctAlnAx < <p{Xx) 
leading to 

(25) ax < -(AlnA)"V(A5), V/j" 1 < A < 1. 

It follows that a is bounded above by the expression in (120]) . To prove that conversely, (TT9]) 
holds when a is given by (|20|) . we use an invariance argument as in the proof of Propo- 
sition E3J Denote: 

X := # x K, 
V :=U x (1,/x), 

51 := {(x, A) G # : x < Ax}, 

5 2 := {(x, A) G # : x < Ax}, 
and define the differential equation: 

<»> (i)=*(^)-( -SlnA 

Then (1191) holds if flD, $) and (S2 C\T>, $) are invariant systems. Given the convexity 
of Si, we can directly apply Theorem 12 .21 to prove the invariance of the system (SiDT), $). 
The tangent cone of <Si at point (x, A) G Si is given by: 

T Sl (x, A) = {(2,77) : (q,r}X - z) ^ 0,Vg G C*, (g, Ax - x) = 0}. 

For any q £ C* such that (g, Ax — x) = 0, by the order-preserving assumption and 
Proposition 13.31 

(q,<j)(Xx)) > (g,0(x)). 
Now, using the expression of a in (1201) . 



(<7, — aA In Ax — <p(x)) ^ (g, —aX In Ax — 0(Ax)) ^ 0. 
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This shows that 

$(x,\) e T Sl (x,X), \/(x,X) eVnSi, 

whence the invariance of (S\ D T>, $). For the invariance of system (S 2 D X>, $), we define 
a bijection on D: 

F(x, A) = (Ax — x, A) 

whose inverse is: 

G(y,K) = {k~ 1 {x + 

If (x(-),A(-)) G X> follows the dynamics of fl2SD, then (y(-),«(0) = F(x(-),A(-)) is the 
solution of the following differential equation: 

(27) f = K ) = ( ~ a M ^ + y ^ + + 

\ Ki J \ Ct/^ 111 K> 

Thus the invariance of system (F(T>) fl F(S 2 ), implies the invariance of system (X> D 
S 2 , $). Note that F(«S 2 ) =Cxl. Therefore by TheoremOthe system (F(D)nF(S 2 ), *) 
is invariant if 

ttfe, «) e 7> (S2) (y, k), V(y, k) g n F(S 2 ). 
The tangent cone of F(S 2 ) at point (?/, k) G F(S 2 ) is given by: 

Tf^jfo, «) = {*: ^0 ) V 9 Gr ) (g,2/}=0}xl. 
Again by the order-preserving assumption, for any q E C* such that (q, y) = 0, 

(q,<KK-\x + y))) > {q,<P(^\x))) 
Using again the expression of a in ( 120]) . 

(g, ^(k, -1 (x))) ^ (g, (aln/«)x) 

because k G Therefore 

(g, -(a In k)(x + y) + + y))) ^ 0, 

which implies 

Hv, K ) e rp (<S2) (y, «), v( y , k) g f(x>) n f(S 2 ), 

whence the invariance of (F(S 2 ) n F(I>), and that of (S 2 nD,$). 

Finally, if a ^ 0, then the set W is invariant (by ffT9]) ). Thus tu{x) = +00 for all 
x eU. Since the closure of U is in the interior of the cone, we conclude that 

the relation f[T§j) holds as well for x G [/J,~ l x, fix]. □ 

3.3. The discrete time case. For completeness, we give in this section the results 
analogous to Proposition 13.31 and Theorem 13.51 for discrete operators, which are of a 
simpler character. In this section we consider a differentiable map F : Cq — » Cq. The 
first proposition characterizes order-preserving maps, its elementary proof is left to the 
reader. 

Proposition 3.8. Let U d Cq be any open convex set. Then F is order-preserving on U 
if and only if 

DF(P)-z^o, \/PeU,ZeC 

Let Q C Cq. The Lipschitz constant of F on Q, denoted by Lip(F; Q), is defined as: 

fort t • (j? r \ rf T (F(P 1 ),F(P 2 )) 

(28) Lrp(F;£):= sup . 

PuPieS a T \f x ,f 2 ) 
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Proposition 3.9. Let Q C Cq be a set such that tQ C Q for any t ^ 1. If F is order- 
preserving on Q , then 

Lip(P; Q) = inf{a : DF(P) ■ P <: aF(P), VP G 6?}. 

Proof. It suffices to prove the equivalence between the following two conditions: 

(a) d T (F(P 1 ),F(P 2 )) ^ ad T (Pi,P2), VP 1; P 2 G Q 

(b) PP(P) • P sC aP(P), VP G £ 

As was pointed out in Remark 1.9 |Nus94] . if F is order-preserving, then Condition (a) 
is true if and only if: 

\- a F(\P) ^ P(P), VP G g, A > 1. 

Condition (b) is a necessary condition (differentiate the above inequality at A = 1). For 
the sufficiency, note that the derivative of the left-hand side is 

\- a -\DF{\P) ■ (AP) - aP(AP)) 

which is always negative semi-definite given that Condition (b) is true. □ 

Remark 3.10. Nussbaum treated the discrete case in |Nus94j . as an intermediate step 
before considering differential equations. Corollary 1.3 there shows that for any open 
subset Q C Co such that for all u, v G Q there exists a piecewise ^ minimal geodesic 
contained in Q (geodesic convexity assumption), the Lipschitz constant of the map F on 
Q satisfies : 

(29) Lip(P; Q) = inf{a : -aF(P) ^ DF(P) ■ Z <: aF(P), VP G £, -P Z ^ P} 

Thus, when the map F is order-preserving, a variant of Proposition 13.91 in which the 
domain Q satisfies the previous geodesic convexity assumption can be easily obtained as 
a corollary of this result. 

4. First applications and illustrations 

In this section, we show that several known contraction results, which were originally 
obtained in |LW94j and [LL07J by means of symplectic semigroups, as well as new ones 
concerning the standard Riccati equation with indefinite coefficients, can be obtained 
readily from Theorem 13.51 The extension of these results to the generalized Riccati 
equation will be dealt with in Section 

4.1. Contraction rate of order-preserving flows on the standard positive cone. 

Let us consider the standard cone C := {x G lR n : X{ ^ 0, 1 ^ i ^ n} in X = W n and 
an order-preserving flow M.'(-) associated to a differentiable function <p : J x C — > M n . 
For a subset U C Co, define the best contraction rate on IA to be the greatest value of a 
satisfying: 

(30) 

driMKx^,^^)) ^ e~ ait ~ s) d T (x 1} x 2 ),Vx u x 2 eW,0^^< t u (s, x x ) A t u (s, x 2 ) 
A direct application of Theorem 13.51 is the following: 

Corollary 4.1 (Compare with |Nus94[ Th. 3.10]). LetU C Cq be an open set satisfying 
\U C U for all A G (0, 1] . For s G J and x dlA define gi(s, x) by: 

(31) 9i(s,x) = -xr 1 [^2^(s,x)x j - (j)i(s,x)\ 

3=1 3 
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then the best contraction rate on U defined in ( 13 Op is given by: 

(32) a = inf{<7j(s, x) : 1 ^ i ^ n, x G U, s G J} 

This should be compared with a result of Nussbaum |Nus94] . who showed that a 
modification of this formula, with an absolute value enclosing each term ^(s,x) for 
i 7^ j, holds for a non-necessarily order-preserving flow. Nussbaum's approach uses 
the fact that the Thompson metric originates from a Finsler structure to determine the 
local contraction rate. However, this leads to different assumptions (see Assumption 
H3.1 in [Nus94j ). In particular, as in the discrete case (see Remark l3.10p . the method 
of [Nus94j requires some form of geodesic convexity assumption, which can be dispensed 
with if the flow is assumed to be order-preserving. For instance, only the special case 
of Corollary 14.11 in which the domain U is geodesically convex can be recovered by the 
method of [Nus94j . 

4.2. Standard Riccati operator. One major application of the above analysis is the 
Riccati operator, arising from the Linear Quadratic (LQ) control problem. Let E be a real 
or complex Hilbert space with inner product (-,•). The set of bounded linear operators 
on E is denoted by End(E). For A G End(E), let A' denote the adjoint of A. The set of 
symmetric bounded linear operators is denoted by Sym(E). A symmetric bounded linear 
operator A is positive semi-definite if (x,Ax) ^ for all x G E. Following [LL07], let 
P(resp. Vq) be the set of positive semi-defmite(resp. positive semi-definite invertible) 
bounded symmetric linear operators of E. Then V is a convex closed pointed cone with 
interior Po(Lemma 9.2 and Proposition 9.5 |LL06j ) and induces the Loewner order 
on Sym(E): 

p <^q <^q- p eV. 

Then we may define the Thompson metric on Vo- This is of course a special case of 
the definition in Section 12.11 Note that equipped with the operator norm, the cone V is 
normal. Therefore the metric space (Vo; dr) is complete(Lemma 5.1 [LL07] ). 
Consider the Riccati differential equation defined on Sym(E): 

(33) P(t) = 4>(t, P) := A(t)'P(t) + P(t)A(t) + D(t) - P(t)E(t)P{t). 

where A : R — > End(E), D : R — > Sym(E), E : R — > Sym(E) are assumed to be continuous 
and bounded applications. The flow associated to fl3"5|) is naturally order-preserving on 
Vo by considering the LQ control problem. One may also verify it using Proposition 13.31 
The least contraction rate of the flow on Vq is the best constant a such that for all 
Pi,P2 e Vo and s ^ 0, 

(34) d T (M*(Pi),M*(P 2 )) ^ e- a{t - s) d T (P u P 2 ), \/s^t< t-p (s, P±) A t-p (s, P 2 ). 
An immediate consequence of Theorem 13.61 is: 

Theorem 4.2. The least contraction rate defined as in (|34|) satisfies: 

(35) a = sup{/3 G R : PH{t)P + D(t) > f3P, Vt ^ 0, P G Vq] 

Remark 4.3. Even if in the statement of Theorem 14.21 we do not require E and D to be 
positive semi-definite, the set of the supremum of which is taken in fl35]) is easily seen to 
be empty as soon as E or D are not positive semi-definite. Hence, the finiteness of the 
constant a in Theorem 14.21 does require E and D to be positive semi-definite and then 
we must have a ^ 0. This shows a dichotomy: either the flow is non-expansive, or it is 
not uniformly Lipschitz. 
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Corollary 4.4 (Theorem 8.5 |LL07j ). We suppose that £>(*),£(*) G V for all t ^ 0. 

Then the least contraction rate is given by: 

« = 2inf A /m((S(t)V2 jD ( t ) S(t )i/2 )//) 

Proof. The best contraction rate is given by: 

a = sup{/3 ^ : P£(t)P + D(t) ^ 0P, Vt ^ 0, P G V } 
Consider all P = XI, then 

a ^ sup{/3 ^ : A 2 £(t) ^ /3AI - £>(t), Vt ^ 0, A > 0} 
If £(t) G P is not invertible, then m(E/7) = 0. Thus 

a ^ sup{/3 ^ : ^ /3A + m(-D(t)/I), Vt ^ 0, A > 0} = 
Now suppose that £(t) G P , Vt ^ 0. In that case, PE(t)P + D(t) ^ /3P if and only if 

S(*)5PE(*)PE(t)2 + E(t)*D(t)E(t)z - /3£(t)sP£(t)i 
= (E(t)5PE(t)^) 2 - /3E(t)5PE(t)5 + E(t)5D(t)E(t)5 
= (E(t)IPE(t)i - f J) 2 + E(t)lp>(t)E(t)l - f / ^ 

Therefore, 

a = sup{/3 ^ : g < 2^m((£(t) 1 / 2 P( t)E(t)V 2 )/I), Vt ^ 0} 
= 2 inf ^m((S(t)V2D(t)S(t)V2)//). 

□ 

The above theorem was proved by Lawson and Lim in [LL07] . Theorem 8.5, using a 
Birkhoff contraction formula of the fractional transformation on symmetric cones. Their 
approach requires the coefficients E(t) and D{t) to be positive semi-definite. By Re- 
mark 14. 3[ this condition is also necessary to the existence of a global contraction rate. 
However, a local contraction may occur even the coefficients are not positive semi-definite. 

We now consider the Riccati equation with constant coefficients (A,D,T,). The com- 
mon fixed point of the flow M t for all t must satisfy the algebraic Riccati equation(ARE) 
equation: 

A'P + PA + D - PEP = 

If E, D G Vq, by Corollary 14.41 and the completeness of the metric space (P j dr) we know 
that the solution of ARE exists and is unique. We next give sufficient conditions for the 
existence of solutions of ARE even when E is not positive semi-definite. Below is a direct 
consequence of Theorem 13.51 

Corollary 4.5. Let Pq G Vq and a G R. The following are equivalent: 

(a) ForallP 1 ,P 2 G (0,P ), 

d T (M t (P u P 2 )) < e- at d T (P 1 ,P 2 ), Vt < t (0)Po) (Pi) At (Q) p o) (P 2 ). 

(b) For all P G (0,P ), 

D + PEP ^ aP. 

In particular, this corollary allows to prove the local contraction property of the Riccati 
equation (I3"3l when E is not positive definite. Let CA,CD,m D ,c^ G R such that: 

A + A' ^ -2c A I, m D I < D ^ c D I, E ^ -c s J. 

The situation considered in the next corollary is motivated by the analysis of a method 
of reduction of the curse of dimensionality introduced by McEneaney |McE07] . This 
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method applies to a control problem in which one can switch between several linear- 
quadratic models. 

Corollary 4.6. Suppose that ca,cd > 0,m D ,cs > and 



<? A ^ c o c s, c s m D > (c A - \lc A - cdCt.) 



then for any A G [ CA c -d c s) ^ /^.^ t/iere is a ^ (m£> — ceA 2 )/A such that for all 

P 1 ,P 2 e(0,XI] 

dr(M t (Pi), d r (P 2 )) < e- at d T (P h P 2 ), Vt ^ 0. 
In particular, there exists a unique solution P to ARE in (0, XI] and for any P G (0, AT], 

d T (M t (P),P) ^ e- at d T (P,P), Vt ^ 0. 



Proo/. Let any A G [f^MEf^ Jm£\ Since 

0(AJ) = A(A + A') + D - A 2 £ (-2Ac A + c D + A 2 c s )J < 0. 

we deduce that the closed set (0, A J] is invariant by the Riccati flow. It is not difficult to 
show that given Ao G (A, J^§) there is a ^ (mjj — ceAq)/Ao such that 

D + P£P ^ aP, VP G (0, A J). 
Indeed, note that a sufficient condition would be: 

m D I - c E P 2 ^ aP, VP G (0, X I) 

which is equivalent to: 

m D - c s Aq ^ «A . 
By Corollary KM for any P 1; P 2 G (0, A/] C (0, A J), 

dr(M t (Pi), Af t (P 2 )) < e- Q *rf T (P 1 ,P 2 ), Vt ^ 

Since the metric space ((0, XI]; dx) is complete, we deduce that there is a unique fixed 
point P G (0, AJ] and all solutions with initial value in (0, XI] converge exponentially to 
P with rate a. □ 

Another interesting case is when £ G P not invertible. In that case, Corollary 14.41 
tells that the least contraction rate on Vq is 0. However, using Corollary 14.51 we can say 
something more about the asymptotic behavior of the trajectories. 

Corollary 4.7. Suppose that ca > 0, mo > and cs = 0. Then for any X ^ ^a, there is 
a > swc/i t/iat t/ie flow is a-contractive on the set (0, XI). In particular, the existence 
and uniqueness of solution P G Po ^° ARE is insured and for any P G Vq, 

d T (M t (P),P) ^ e- at d T (P,P), Vt ^ 0. 

where a = min(m(J/P), ^)mfl. 



We leave the proof to the reader, which is similar to the one of Corollary 14.61 
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5. Application to stochastic Riccati differential equations 



In the sequel, denote by S n the set of n-dimensional real symmetric matrices. Observe 
that equipped with the canonical inner product 

(A,B) = Tr(AB), V7L, B G S n , 

S n is a real Hilbert space. The subset of all positive semi-definite matrices S+ C S ra forms 
a closed pointed convex cone. All the above results can then be applied here considering 
S™ as C and the set of all positive definite matrices S™ + as C . Note that here C* = C. We 
shall use the notation ^ (and 3>) for the (strict) Loewner order, and dx for the Thompson 
metric induced by S™ (see Section 12.11) . 

5.1. Stochastic LQ problem and GRDE. Consider the following stochastic linear 
quadratic optimal control problem: 

v(s, y) miii E f \x(t)' Q(t)x(t) + 2u'(t)L(t)x(t) + u(t)' R(t)u(t)]dt + E[x(T)'Gx(T)) 
«(*) Js 

dx(t) = (A(t)x(t) + B(t)u(t))dt + (C(t)x(t) + D(t)u(t))dW(t), Vt G [s, T], 
x(s) = y. 



s.t. 



where the functions appearing above satisfy: 

A(-),C(-) G L°° nC°(0,T;M nxn ), B(-), D(-), L(-) G L°° fl C°(0, T; M nxfc ), 
Q(-) G L°° PI C°(0, T; S"), #(•) G L°° n C°(0, T; S fe ). 

Here is a standard Brownian motion defined on a complete probability space. We 
refer the reader to [YZ99J Chapter 6, for the precise definition of this control problem. 
In [YZ99] . the above functions are only assumed to be bounded. In our case, the conti- 
nuity is necessary to apply the previous results. 

The above stochastic LQ control problem over the time interval [s,T] is solvable, i.e., 
admits an optimal control for all y G M. n if the solution of the following constrained 
differential matrix equation exists: 



(36) 



P + A'P + PA + C'PC + Q = 

(PB + C'PD + L')(R + D'PD)~ 1 (B'P + D'PC + L), t G [s,T] 
P{T) = G 

R(t) + D(t)'P{t)D{t)^>0, te[s,T\ 



which we refer to as generalized Riccati differential equation (GRDE). In that case, the 
value function of the optimal control problem is given by 

(37) v(s,y) = y'P(s)y. 

5.2. GRDE with semi-definite weighting matrices. The solvability of the GRDE (|36|) 
with indefinite matrix coefficients has been treated by Chen, Moore, Rami, and Zhou 
in [RCMZOfb] . In order to apply our previous results, we only consider the case: 



(38) ( ^ ^ \> 0, kerR(t)nkerD(t) = {0}, Vt G [0, 



n 



so that the function 



, v 4>(t,P)= PA + A'P + C'PC + Q- 

[6J) (B'P + D'PC + L)'(R + D f PD)-\B f P + D'PC + L) 
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is well defined on [0, +00) x S™ + and satisfies the assumptions made at the beginning of 
section [3J We are going to apply the preceding results to show the monotonicity and the 
non-expansiveness of the GRDE differential equation defined on S™ + : 



(40) 



P = (f>(t,P), 
P(0) = G 



Proposition 5.1. Assume that fl38|) holds. Then the flow associated to [JCty is order- 
preserving and non- expansive on S™ + . 

Proposition 15.11 could be derived by exploiting the relation between the solution of 
the Riccati equation and the value function of the stochastic control problem (see (1571) ). 
Here we choose to prove it from the infinitesimal characterizations of Proposition 13.31 and 
Theorem 13.61 



Proof. By Proposition I3.3[ if suffices to prove that for any P E S™ + , any Q, Z E S" such 
that (Q, Z) = 0: 

(Q,D<f> t (P)Z) >0. 

Indeed, 

D<f) t (P)Z = ZA(t) + A{t)'Z + C(t)'ZC(t) - (B(t)'Z + D{t)'ZC{t))'N t {P) 
-N t {P)'{B{t)'Z + D{t)'ZC{t)) + N t (P)'D(t)'ZD(t)N t (P) 

where N t (P) = (R(t) + D(t)' 'PD^))' 1 ^^) 1 P + D{t)'PC{t) + L{t)). Remark that if 
Q, Z E S'l and (Q, Z) = then QZ = 0. Therefore, 

(Q, D<t> t {P)Z) = (Q, C{t)'ZC{t) - C(t)'ZD(t)N t (P) - N t (P)'D(t)'ZC(t) 
+N t (P)'D(tyZD(t)N t (P)} 
= (Q, (C(t) - D{t)N t {P))'z(C{t) - D(t)N t (P))) > 0. 

Now for non-expansiveness, by Theorem 13.61 it remains to verify that for any P E S™ + 
and any t E [0, T], 

D<p t (P)P-<f)(t,P)^0. 



Indeed, 
(41) 



D<p t (P)P-<P(t,P) 

= -Q(t) + N t (P)'L(t) + L'(t)N t (P) - N t (P)'R(t)N t (P) 



where H t (P)' = ( / -N t (P)' ). □ 

Remark 5.2. A fundamental discrepancy with the standard Riccati equation is that the 
flow of the generalized Riccati equation is not a global contraction. This is because that 
there is no a > such that the condition 



D(f> t (P)P - <f)(t, P) ^ -aP, VP E S r , 



■+' 



which by Theorem 13.61 is necessary to the global contraction property of the flow, is 
satisfied. However, we shall see in the next section that a local contraction property does 
hold. 
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5.3. Asymptotic behavior of GRDE. We are going to investigate the behavior of 
the GRDE flow as time horizon goes to infinity. All the matrices A, B, C, D, L, Q, R are 
assumed to be constant. First we show a local contraction property under the condition 

Q V 



(42) y L R j»0. 

More precisely, 

Theorem 5.3. Assume that ( 142]) holds. LetU C S™ + be an open set such that \U C U 
for all A G (0,1]. Assume that there is P$ G S™ + such that U C (0, Pq] an d let a = 
m(Q - L'R^L/Po), then for all P u P 2 eU, 

dr(M(Pi), M t (P 2 )) ^ e- at d T (Pi, P 2 ), 0^t< t u {Pi) A * M (P a ) 

Proof. By applying Theorem 13.51 we need to prove 

D(j){P)P - <j>{P) < -aP, VP G U 

Indeed, for all P G U, 

Q-aP- L'R- X L ^Q-aP - L'R~ l L > 0. 
Besides, the previous calculus yields 

P0(P)P - 0(P) + aP 

(43) =p(py(:^+- p :^')^) 

where H(P)' = ( / -N(P)' ) and JV(P) = (P + D'PD)- 1 (B'P + P/PC + L). By 
Schur's complement lemma, we get 

P0(P)P - 0(P) ^ -aP, VP G W. 

□ 

The fixed point of the GRDE flow associated to (HOI) , if it exists, satisfies the so-called 
general algebraic Riccati equation (GARE): 



(44) 



0(P) = 0. 

P + D'PD > 



where 0(P) := A'P+PA+C'PC+Q-(B'P+D'PC+L)'(R+D'PDy 1 (B'P+D'PC+L). 
The existence of solutions of GARE and the asymptotic behavior of the GRDE flow have 
been studied in [RCMZOlb] and [RZ00J. The authors assumed the following mean-square 
stabilizability condition: 

Definition 5.4 (Definition 4.1 [RZOO] ). The system of matrices (A, B,C, D) is said to 
be mean-square stabilizable if there exists a control law of feedback form 

u{t) = Kx(t), 

where K is a constant matrix, such that for every initial (to,xo), the closed loop system 

dx(t) = (A + BK)x(t)dt + (C + DK)x{t)dW{t) 
x(0) = xo 



satisfies 



lim E[x(t)'x(t)} = 

t— H-oo 
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Under the mean-square stabilizability assumption, they established a necessary and 
sufficient condition for the existence of a solution. To make a comparison, let us first 
quote their theorem: 

Theorem 5.5 (Theorem 4.1 [RCMZOlb]). Under the mean-square stabilizability assump- 
tion, there exists a solution of the GARE ftJ4\ ) if and only if there exists Pq G S n such 
that 

0(P O ) > 0, R + D'P D > 
Moreover, for any such Pq, the solution P(t) of ( f^0| j with initial condition P(0) = Po 
converges to a solution to the GARE as t — )■ oo. 

It follows directly from the above theorem that under the mean-square stabilizability 
assumption, if (j42j) is true, then there must be a solution to the GARE f l44"j) . We next 
show a necessary and sufficient condition for the existence of a stable solution without 
the mean-square stabilizability assumption. 

Theorem 5.6. Assume that the condition ( 142|) holds. Then, the GARE admits a solution 
P G S'1 + if and only if there exists P G S" + such that: 

(45) 0(P O ) ^ 0. 
In that case, for any P G S™ + : 

d T (M t (P),P) ^ e- at d T (P,P), Vt ^ 0, 

where 

1 _ e -d T (P,P) 
01 > d T (PP) m((g " L ' R ~ lL V P ) > °- 
In particular, the solution is unique in S" + . 

Proof. If P G S™ , is a solution of the GARE, then ( 14"5]) is satisfied by considering Po — P ■ 
Conversely, note that if 4>(Pq) ^ for some Po G S+ + , then (0, Po] is an invariant set. 
Consider the open set U = (0, Pq + I). By Theorem I5.3[ there is a > such that for all 
Pi, P 2 G (0,P ] C W, we have: 

d T {M t (P x ), M t (P 2 )) < e- Qi rf T (Pi, P 2 ), V0 ^ t < t u (Pi) A t w (P 2 ), 

Since [0, Po] C U is invariant, we have that tw(Pi), t^(P 2 ) = +oo. Thus the flow M t is 
contractive in the complete metric space ((0, P ], d T ). There must be a unique fixed point 
P G (0, Po] such that <f>(P) = 0. Next, assuming the existence of a solution P G Co to the 
GARE, we apply Theorem 13 .71 to obtain the rate of convergence. A basic calculus yields: 

A-V(AP) 

= (B'P + D'PC)'((R + D'PD)- 1 - {X^R + D'PD)- X ){B'P + D'PC) + (A -1 - 1)Q 
where C = C- DR^L and Q = Q - L'Rr l L. Therefore, if A > 1, then 

(46) A-V(AP) ^ (A" 1 - 1)Q 
and 

(47) A0(A- 1 P) ^ (A - 1)Q. 

Now for any P G S™ + ^ P, let // = e^ p '^ and a = ±^m(Q/P) > 0. Then 

(A" 1 - 1)Q ^ -a(lnA)P, (A - 1)Q > a(lnA)P, VA G (1,/i) 
and (g6]) and (07} lead to: 

aln(A)P < A0(A~ 1 P), cdn(A)P ^ -A~V(AP), VAg(1,/x). 
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Thus, 

(48) 0<a^ inf m((-(A In A)~V(AP))/P) • 

At~ 1 <A</x 

By virtue of (14"8"]1 and Theorem 13.71 we have 

d T (M t (P),P) < e- a 'rf T (P,P), Vt ^ 0. 

□ 

5.4. Discrete Generalized Riccati operator. The linear quadratic stochastic control 
problem has a discrete time analogue |RCZ01j . which leads to the generalized discrete 
Riccati operator P : S n -> S n : 

(49) 

P(P) = A PA + C'PC + Q - (B'PA + D'PC)'(R + P'PP + D'PD)-\B'PA + P/PC) 

where A,C E W nxn , B,D G R nxm and Q,R e S n . We assume that Q > and R > 0. 
Then by applying the Schur complement condition for positive defmiteness, one can prove 
that F sends S" + to itself. Note that when C = D = 0, we recover the standard Riccati 
operator: 

(50) T(P) = A' PA + Q- A'PB(R + B'PB)- 1 B'PA. 

The object of this section is to get the Lipschitz constant of F on S™ + (see (128]) ). First 
we show that this operator is order-preserving on , . 

Proposition 5.7. The operator F is order-preserving on S? + . 

Proof. Let any P G S™ + and Z 6 S". A simple calculus show that: 

DF(P) ■ Z = (A- BN)'Z{A - BN) + (C - DN)'Z(C - DN) ^ 

where N = (R + B'PB + D'PD)-\B'PA + D'PC). By Proposition P is order- 
preserving on S+ + . □ 



Next we apply Proposition 13.91 to get: 
(51) Lip(P; S^ + ) = inf {a ^ : DF(P) ■ P <: aP(P), VP G S™ + }. 

The following two lemmas will be useful. 

Lemma 5.8. Let = (^)w fre a rank factorization (so that the last two factors 
have maximal column and row rank, respectively). Then the operator F defined in Qjty 
satisfies: 

(52) 

P(P) = A' PA + C'PC + Q- (B'PA + D'PC)'(R + B'PB + D'PD)-\B'PA + D'PC) 
where R = (WR^W)- 1 . 

Proof. To simplify the notation, denote X(P) = BPB + D'PD. Notice that since the 
matrix is of full column rank, X(P) is invertible for all P G S™ + . It follows from (T4"9"j) 
that: 

P(P) = A' PA + C'PC + Q- (B'PA + D'PC)'W(R + W'X(P)W)' l W'(B'PA + D'PC) 

Now appealing to the Woodbury matrix identity, we obtain: 

(53) 

W(R + X(P)W)~ 1 W = W(R~ l - R- 1 W'(X(Py l + WR^W'^WR-^W 

= WR~ l W - WR~ l W'(X(P)- 1 + WR^W'^WR^W 
= ((WR^W)- 1 + X(P))- 1 
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from which we get (I52p . □ 
Lemma 5.9. Let 5^2, then 

(54) X - X(R + Xy\5R + X)(R + X)~ X X < — — - , VI G S n + 

Proof. Let any X e S" . Since X commutes with /, we have that: 

X -X(I + X)- 1 (5I + X)(I + X)- X X 

= (/ + X)-\X(I + X) 2 - X\5I + X)){I + X)- 1 

= (/ + X)-\{2 - 5)X 2 + X- ^(1 + Xf){I + X)- 1 + ^1 

= -(I + XY\{25 - 3)X - If {I + X)" 1 + 

^ 4(5-1) J ' 

To obtain (|54p . it suffices to notice that: 

Rr*(X - X(R + X)-\5R + X)(R + X)~ 1 X)R~^ 
= Y-Y(I + Y)- 1 {5I + Y)(I + K) -1 Y 

where Y = iHxiH. □ 
Proposition 5.10. The operator F is non- expansive: Lip(F; S™ + ) ^ 1. Let 



B )= B 
D \ D 



W 



be a rank factorization. Then a necessary and sufficient condition to have Lip(F; S™ + ) < 1 
is that there is a matrix S such that: 

(55, ( A c)-( S B )S. 



In that case, 



++ (1 + ^l + M(S'RS/Q)f 



1WW 



where R = (WR^W 

Proof. Lemma 15781 implies that it is sufficient to prove the proposition for the case W = I, 
i.e. when (^) is of full column rank. A simple calculus shows that: 

DF(P) ■ P - aF(P) = (1 - a) (A' PA + C'PC) - aQ 

-(1 - a)N{P)'{R + X(P))- 1 X(P) 
-N(P)'(R + X(P)y l R(R + xIp))- 1 N(P) 

where 

N(P) = B'PA + D'PC, X(P) = B'PB + D'PD. 

Then it is evident that Lip(F; S" + ) < 1. Now let S G R nxm such that §B) holds. Then 
N(P) = X(P)S, A' PA + C'PC = S'X(P)S and 

DF(P)-P-aF(P) = (l-a)S'X(P)S - aQ 

-(1 - a)S'X(P)(R + X(P))- l X(P)S 
-S'X{P)'{R + X{P))- l R{R + X{P))- l X{P)S 

To simplify the notation, let X := X(P) and 5 := |z^, then 

DF{P) ■ P - aF(P) = (1 - a)S'(X - X(R + Xy^dR + X)(R + X)~ 1 X)S - aQ. 

22 



By Lemma [5.91 

X - X(R + X)- 1 ^ + X)(R + X)- l X <: = ^R, VX G S™ . 

Therefore 

(1 - a)S'(X - X(R + X)-\5R + X)(R + X)~ 1 X)S ^ ^~ a ^ S'RS. 

Consequently if a is such that: tjz^W = M(S'RS/Q), then 

DF(P) ■ P - aF(P) < 0, VP G S™ . 
Together with ()5ip this shows that 

M{S'RS/Q) 



Lip(P; S" ) < 



i + v/T+ M(S'i?s/g)) 2 



Next we prove the necessity of condition fl55l) . Remember that since the matrix (^) has 
full rank, X(P) is always invertible for P G SIL_. Besides, there is a < 1 such that 

DF(P) ■ P - aF(P) < 0, VP G S r | + 
if and only if for any P G S™ + , 
(ATM + C'PC) - -^Q - N(P)'(R + X(P))- 1 (f5f R + X(P))(P + X(P))" 1 X(P) ^ 0. 
That is, for any P G S" + and A > 0, 
(A'PA + C'PC) - f^Q - X(P)'(±P + X(P))-\^R + X(P))(±P + X(P))^N(P) 
Letting A go to infinity, by continuity, we obtain that: 

(A'PA + C'PC) - N{P)'X{P))- l N{P) <: 0. 
The above expression is the Schur complement of the positive semi-definite matrix 

(B'PB + D'PD B'PA + D'PC\_(B'\ { , ( D> \ , , 

\ A'PB + C'PD A'PA + C'PC ) ~ \ A' J v > + \^ C ) v u ^ )• 

Therefore for any x G W 1 there is u G IR m such that 

,( u\ (B'PB + D'PD B'PA + D'PC \ ( uV_ 
W ^ x y ' ^ A'PP + C'PD A'PA + C'PC )\x ) ' ~ 

That is, for any iGl™ there is u G W 11 such that: 

B A \ f u 
D C )\x 

This is equivalent to say that there is S G ]R mx " such that: 

A \ = ( BS 
C ) \DS / 

□ 

The contraction rate of the standard discrete Riccati operator T : S™ + — > S™ + can now 
be recovered corollary: 
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Corollary 5.11 (Compare with [LL08J). The standard Riccati operator T defined in f)50p 
is non-expansive: Lip(T; S™ + ) ^ 1. A necessary and sufficient condition to have the strict 
contraction property is that the matrix B is of full row rank. In that case, let B = BW 
be a rank factorization, then 

L 1P (T;S« ++ K M{S ' RS/Q) < 1 

++ (1 + ^1 + M(S'RS/Q)y 

where S = B~ X A and R = (WR^W')' 1 . 

Remark 5.12. Condition (1551) leads to a formal argument explaining why strict global con- 
traction cannot be hoped for the GRDE flow. Indeed, we can approximate the continuous- 
time LQ control problem in Section [5] over a small time horizon e by the following one-step 
discrete-time stochastic linear quadratic control problem: 

min E((x , eQx ) + (u, eRu) + (x t , Gx e )) 

u£l m 

s.t. x e — (I + eA)xo + eBu + (^Cx + y/eDu)w 

where w ~ A/"(0, 1). Without loss of generality, we suppose that (^) is of full column 
rank. If a strict contraction result was valid for the continuous time system, we would 
expect the same to be true for its discrete approximation if e is sufficiently small. However, 
the strict global contraction condition requires the existence of 5* such that: 

I + eA \ _ f eBS 
yfeC )~\ yfeDS 

which can not hold for a set of e converging to if C and D are not zero. 

6. LOSS OF NON-EXPANSIVENESS OF THE GRDE FLOW IN OTHER INVARIANT 

FlNSLER METRICS 

The standard Riccati flow is known to be a contraction in the standard Riemannian 
metric |Bou93] , and more generally in any invariant Finsler metric (with the same bound 
on the contraction rate) [LL08]. We next construct an explicit counter example showing 
that the Thompson metric is essentially the only invariant Finsler metric in which the 
GRDE Riccati flow is non-expansive. 

6.1. Preliminary results. We first recall the definition of symmetric gauge functions 
and of the associated invariant Finsler metrics on the interior of the cone of positive 
definite matrices. Then, we will show some conditions that are necessary for an order- 
preserving flow to be non-expansive in a given metric of this kind. 

Definition 6.1 (Symmetric gauge function). A symmetric gauge function v : W 1 — > K is 
a convex, positively homogeneous of degree 1 function such that for any permutation a, 

u(X 1 , ■ ■ ■ , A n ) = • • • , \K(n)\), V A = (Ai, • • • , A n ) E R n . 

The next lemma collects several useful properties of subdifferentials of symmetric gauge 
function (see |Roc70] for more background on subdifferentials). The straightforward proof 
is left to the reader. 

Lemma 6.2. Let v : W 1 — >■ R be a symmetric gauge function. The following properties 
hold: 

1 For all AGl n and fie du(X), 

fiiXi ^0, Vi = 1, ■ • • ,n. 
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2 For all A G K n and fie dv{X), 

(fi,X)=v(X). 

3 For all A, A' G R n and fi G du(X),fi' G dv(X'), 

(fi — fi', A) ^ 0. 

For every symmetric gauge function u, we define a spectral function z> : S n — > M: 

*(P) = u(X(P)). 
where A(P) is the vector of eigenvalues of P. 

Theorem 6.3 (Lcw96j). If v is a symmetric gauge function, then v is a convex function 
on S n . Moreover, Z G d£>(P) if and only if there exists y G du(X(P)) such that: 

Z = Vdizg(y)V T , 

where V is the unitary matrix such that P = V dia.g(X(P))V T . 

Following |Bha03| . |LL08] and |ACS 00j. we define a metric on S™ + as follows, 

d u (P,Q) = 0(\og(p- 1 / 2 QP- 1 / 2 )). 

It coincides with the Finsler metric obtained by thinking of S™ , as a manifold and taking 

\\dQ\\ P = u(p-$(dQ)p-t) 

as the length of an infinitesimal displacement in the tangent space at point P. This 
metric is invariant by the canonical action on the linear group on S" + . 
We shall consider specially, as in [ACS00J, the p-norm function: 



<A) = ||A|| p = (^|A,n 1 



/p 



n , 

i=i 

so that the metric d v is the Thompson metric for p = +oo and the Riemannian metric 
for p = 2. 

Lemma 6.4. Let v be a symmetric gauge function and d v be the associated metric on 
S™ + . Let M : S™ — > S™ be a differentiable function such that: 

(57) d u (M(P), M(Q)) <: d u (P, Q), V P, Q G S^ + , 

then 

v(M(Py ll2 (DM(P) ■ Z)M(P)- 1/2 ) < U(P- 1/2 ZP- 1/2 ), V P G S™ + , Z G S" . 

Proof. Let any P G S™ + and Z G S n . There exists S > such that for any ^ e ^ 5, 
P + eZ G S" + . By <$57§ and the definition of d„: 

u\og(M(Py l / 2 M(P + eZ)M(P)- 1 / 2 ) ^ u\og(p- l / 2 (P + eZ)p- 1 / 2 ). 

Divide the two sides by e and take the limit: 

, z>log( M( PY 112 [ M( P + eZ)M(P)- 1 / 2 ) , v \og(P~ l l 2 (P + eZ)P~ 1 / 2 ) 
lim lim 

In view of homogeneity and continuity of the function z>, 

„ \og(M(Py 1 / 2 M(P + eZ)M(P)- 1 / 2 ), „ logP-^fP + eZlP- 1 ^ 
z/(lim ) < z/(lim ) 
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The matrix function log is different iable at /: 

(58) Km k * I tY?- U =0. 

\\u\\^o \\U\\ 

Hence by chain rule: 

P(M(P)- l ' 2 (DM(P) ■ Z)M(P)- 1 ' 2 ) ^ u{p-^ 2 ZP-^ 2 ) 



□ 



(59) 



We consider the following time independent differential equation: 

x(t) = $(x(t)), 



x(s) = x . 



where $ is differentiable on S+ + . We assume that the associated flow M.(-) : (0, +oo) x 
S" + — > S n leaves S" + invariant and is globally defined. 

Lemma 6.5. Let v be a symmetric gauge function. If there exists e > such that for 
any ^ t ^ t, 

(60) i>(Mt(/)- 1/2 (D M t {I) ■ Z)M t {I)- 1 / 2 ) < v(Z), V Z e S n . 

then 

(diag(/i), D $(/) • diag(A) - diag(A)$(/)) < 0, VA G W 1 , \i G dv(X). 
Proof. Let any Z G S". For readability, denote 

P t :=M t {I), H t :=Pl'\ Q t := P~ 1/2 , G t := Pf\ 

and 

U t :=BM t (I)-Z, J t :=U t G t , K t = Q t J t H t . 
The derivative of J t with respect to t is: 

/^-|\ Jt = UtGt — UtGtPtGt 

{ } = (D $(P t ) • l/ t )Gt - U t G t $(Pt)G t . 

The derivative of ii't with respect to t is: 

tf t =Q t j t H t + Q t JtH t + Q t J t H t 

= Q t J t H t — Q t H t Q t J t H t + Q t J t H t 
= Q t J t H t - Q t H t K t + K t Q t H t 

Hence, 

j t | t=0 = D$(/)-Z-Z$(/), 

and 

(62) K t \ t=0 = D $(/) ■ Z - - (4| t=0 )£ + Z(#*|t=o)- 

By Theorem 16.31 the right derivative of the function u(K t ) with respect to t exists: 
0(K t )' + = sup (y,K t ) = sup (Vdmg(v)V',k t ). 

yedv(K(t)) nedv(X),VV'=I 
V'K t V=d\a.g{\) 

Since 

u(K t )^u(K ), te[0,6), 
the right derivative at t = must be negative: 

*(#t)+|t=o < 0. 
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Namely, 

sup (Vdiag(//)F',D$(I) • Z - Z$(J) - (H t \t=o)Z + Z(H t \ t=0 )) *C 0. 

fi£dv(\),VV'=I 
V ZF=diag(A) 

Note that for any unitary matrix V such that V'ZV = diag(A), we have 

(Vdiag(fi)V',(H t \ t=0 )Z) = (Vdiag(fi)V',Z(H t \ t=0 )) 

Hence by taking Z = diag(A) and V = I, we obtain a necessary condition of (1601) : 

(diag(/i), D $(/) ■ diag(A) - diag(A)$(/)) < 

for all A G R n and fi G dv(X). □ 

The above two lemmas lead to the following conclusion: 

Proposition 6.6. If the flow M.(-) : (0, +oo) x S™ + — > S™ + is non- expansive in the 
metric d v , then, 

(63) (diag(/i), D $(/) ■ diag(A) - diag(A)$(/)) < 

for all A G W 1 and \i G dv(X). 

6.2. The counter example. We finally arrive at the announced counter example: we 
give a system of parameters (A, B, C, D, L, Q, R) such that the corresponding $ of GRDE 
does not satisfy the necessary condition (16"3"1) of non-expansiveness in any Finsler metric 
other than the Thompson metric. 
Recall that 

$(P) = A'P + PA + C'PC + Q- (B'P + D'PC + L)'(R + D'PD)- 1 {B'P + D'PC + L). 

Let I n denote the n-dimensional identity matrix and e = (ei, • • • , e„_i)' G IR n_1 be a 
vector. The parameters are chosen as follows: 



A_ T R _ / (e - VI - e)4-i \ _ / (1 + VI - e)/ n _i e 

4 - * - 1 -(vr^iy e J' G -U vr 



and 



D = (VT^e)I n , L = 0, R = el n , Q = el n 

to make 

R + D'D = I n , B' + D'C = I n , C-D = ( I Q n - 1 6 Q ) 

An elementary calculus yields 
D${I) ■ Z - Z${I) 

= A'Z + ZA + C'ZC - (B'Z + D'ZC)'(R + D'D)- 1 (B' + D'C) 
-(B 1 + D'C)(R + D'D)- 1 (B'Z + D'ZC) 
+ (B' + £>'C)(i2 + D'R>y 1 D l ZD{R + D'D)-\B' + D'C) 
-Z(A' + A + C'C + Q - (5' + D'C)'{R + D'D)-\B' + D'C)) 

= 2Z + C'ZC - (B'Z + D'ZC)' - [B'Z + D'ZC) + D'ZD - Z - ZC'C - Z'Q 

= Z + {C - D)'Z{C — D) — B'Z - ZB' - ZC'C - ZQ 

Now let any A = (Ai, • • ■ , A n ) G W l and /i G dv(\). Then 

(diag(//),D$(/) ■ diag(A) - diag(A)$(/)> 
= -2e(//, A) + /x„(-A n |e| 2 + V?) 

Recall that 

(fi, A) = z/(A), V//G3i/(A). 
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So if there is any A G R n and /x G <9i/(A) such that 

n-l 

f^n{-K\e\ 2 + ^X i e 2 l ) > 0, 
i=i 

then there always exists e G (0, 1) such that 

(diag(/x),D$(I) • diag(A) - diag(A)$(/)) > 0. 
Finally we need a lemma to conclude: 

Lemma 6.7. If for all A G M n , /x G cxV(A) and e G M n_1 we /zai>e 

n-l 

M-AJef + ^A^ 2 ) ^ 0, 

i=l 

t/zen 

u(Xi, • • • , A n ) = cmax |Aj| 

i 

for some constant c > 0. 

Proof. First consider e = the z-th standard basis vector of M™ _1 for all % — 1, • • • , n — 1. 
We see that 

A*n(-A n + Ai) < 0, Vz' = l,---,n-l 
for all A = (Ai, ■ • • , A n ) G M n and /x = (/xi, • ■ ■ , /x n ) G dv(X). By the symmetric property 
of z/, this implies actually 

(64) /x^-A,- + Ai) ^ 0, Vz',j = l,--- ,n 

Therefore, for any A ^ if A-,- = then /Xj = for all /x G dv(X). Next, let any 
z G {1, • • ■ , n}, consider the following set 

Ai := {A ^ : Ai = A 2 = • ■ ■ = A, > X i+1 >...Xn> 0}. 

Let any A G Aj and /x G <9z^(A). By Property 1 in Lemma 16.21 /x ^ 0. Using (j64j) . we 
know that: 

/Xj ^ 0, Vj = z + 1, . . . , n. 

Hence, 

Hj = 0, Vj = z + 1, . . . , n. 
Now let any A 1 , A 2 G Aj and /x 1 G ^(A 1 ), xz 2 G <9z/(A 2 ). By Property 3 in Lemma IBT2"l 

(/x 1 - /x 2 , A 1 ) ^ 0. 

It follows that 

i i 

i=i i=i 

We deduce that J^* =1 /xj = X^=iA*f- Hence there is a constant q ^ such that 

3 3 

u(X) = (/x, A) = y~]nj\j = Xi^fMj = CjAi, VA G A i; /x g 0z/(A). 

i=l i=l 

It remains to prove that q = ci for all z = 1, . . . , n. To see this, again we use Property 3 
in Lemma [6.21 First consider A = (1, . . . , 1, 0, . . . , 0) G Aj and any /x G <9^(A), then 

i 

(/x - (ci, 0, . . . , 0)', A) = £ /Xj - ci = Cj - ci ^ 0. 

j'=i 
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On the other hand, for all A 1 G Ai 



(( Cl , o, . . . , oy - fi, a 1 ) = ( Cl - nM - > 



J'=2 



This implies c\ = /ij = q for alH = 1, . . . 
The proof of Theorem 11.31 is now complete. 



n. 



□ 
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